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Superlattice of spin vortices has been proposed in an earlier article as an alternative to the stripe 
interpretation of spin modulations in 1/8-doped lanthanum cuprates. The present article addresses 
several additional characteristics of the spin vortex lattice, namely: (i) the nature of extended charge 
states; (ii) the position of neutron spin peaks as a function of doping; and (iii) the absence of higher 
order spin harmonics. All these characteristics afford straightforward connection to the experimen- 
tal results produced by angle-resolved photoemission spectroscopy, resistivity measurements, and 
neutron scattering. 



The tendency of high-temperature cuprate supercon- 
ductors towards the formation of electronic inhomo- 
geneities is well known and discussed (see, e.g., Ref. 1 
and references therein). Yet the generic outcome of this 
tendency and its relation to the mechanism of supercon- 
ductivity remains unclear. One can only be certain that 
the implications of the inhomogeneous structures would 
be strongly dependent on their dimensionality. A group 
of materials that clearly show inhomogeneous static spin 
response in the bulk is 1/8-doped lanthanum cuprates. 
At present, the majority viewpoint is that this response 
originates from a pattern of one-dimensional spin and 
charge modulations called "stripes"—. This author had 
previously voiced doubt in the stripe interpretation in 
Refsi^. The original consideration in the above refer- 
ences was based on a two-dimensional spin-collinear al- 
ternative to stripes called "grid" . A later experiment 
has, however, produced a result incompatible with the 
grid superstructure 5 , but it did not rule out the possibil- 
ity of a non- collinear two-dimensional superstructure. In 
Ref£, this author has introduced a minimal realization 
of such a superstructure referred to below as "spin vor- 
tex lattice" (see Fig. [IJ , and further argued that it still 
constitutes a viable interpretation of neutron scattering 
experiments in 1/ 8-doped lanthanum cuprates — in par- 
ticular, more straightforwardly than grid reproducing the 
positions of the charge peaks. [Reference 6 - refers to the 
same objects as "magnetic vortices." Here I decided to 
change the terminology to avoid false associations with 
superconducting vortices.] 

Further adding to the difficulties of the stripe inter- 
pretation, a very recent experiment 7 - has reported static 
spin modulation in La1.g7Sro.13Cuo.99Feo.01O4 in the low 
temperature orthorhombic (LTO) phase rather than the 
low-temperature tetragonal (LTT) phase — the latter be- 
ing the cornerstone in the system arguments supporting 
stripes™. The above finding is consistent with the view 
of the author—, that it is not the geometry of the LTT 
phase as such but rather the proximity to the LTO-LTT 
phase transition that stabilizes the static spin response. 

The purpose of this article is to communicate several 
additional facts and considerations related to the spin 
vortex scenario. The issues to be considered are: (i) Ex- 
tended charge states; (ii) Position of neutron spin peaks 



as a function of doping; (iii) Absence of higher order spin 
harmonics. 

The spin vortex lattice shown in Fig. [1] is obtained for 
a square lattice of spins by a coherent superposition of 
two orthogonal spin harmonics with transverse in-plane 
polarizations: 



S« = (-1)*+' 



^ J cos(qi • nj + tpi) 
S ° J cos(q 2 • r y + ip 2 ) 



, (1) 



where qi = (q , 0) and q 2 = (0, go) are the wave vectors 
of the two harmonics with qo = a is the lattice period, 
So is the polarization amplitude of each harmonic, ry are 
the radius- vectors of the lattice sites with the origin cho- 
sen on the one of the sites, and <p%, tf2 are the phases. 
Fig. Q]i corresponds to ipi = ip% = (site-centered vor- 
tices), Fig. [TJd to ip 1 = if2 = \q$a (plaquette-centered 
vortices). Doping level 1/8 implies two doped holes per 
spin vortex. 

This author has not attempted to derive the above 
structure from a microscopic treatment but rather looked 
at the experimental phenomenology. One can, however, 
learn about various theoretical considerations leading to 
the formation of spin vortices from Refs^— . It remains 
a challenge though to understand how two holes would 
organize themselves within a single spin vortex. 

Extended charge states. 

Spin vortex lattice commensurate with the underlying 
atomic lattice exhibits a network of diagonal lines, along 
which spins form perfect spirals — see Fig. Q] These spi- 
ral lines can host gapless extended states of charge carri- 
ers and, therefore, explain the results of the recent angle- 
resolved photoemission spectroscopy (ARPES) study 16 , 
which indicated that Lai.875Bao.i2sCu04 exhibits the 
usual pseudogap structure with gapless direction at 45 
degrees with respect to the principal crystal axes. This 
ARPES result poses a serious challenge to the stripe pic- 
ture. 

The network of spirals can also explain why the resis- 
tivity of the sample drops on cooling below the tempera- 
ture, where the system starts exhibiting static magnetic 
response 1 - 7 . Indeed, the static transition can be inter- 
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FIG. 1: (Color online) Two examples of commensurate spin 
vortex lattices: (a) Vortex cores site-centered, (b) Vortex 
cores plaquette-centered. Straight lines indicate the direc- 
tions, where the structures exhibit perfect spirals. 



preted as the crystallization of vortices, which exist in 
the "liquid" form above the transition temperature. The 
liquid of spin vortices distorts the above network of spi- 
rals and, therefore, exhibits higher resistivity. The resis- 
tivity drop is then explained by a sudden formation of a 
perfect spiral network. 

Neutron peak splitting as a function of doping. 

If one assumes, as normally done in the framework of 
the stripe interpretation, that the static spin superstruc- 
ture observed at 1 /8 doping also exists at other dopings in 
the fluctuating form and thus responsible for the inelastic 
neutron scattering peaks at (-, - ± —5), (£ ± ^-5, £), 
then the spin vortex picture leads to the dependence of 
the dimensionless splitting parameter 8 on the doping 
concentration x (per in-plane Cu atom), which is differ- 
ent from 8 = x obtained for quater-filled stripes. 

Spin vortex lattice contains two holes per vortex. If 
the vortex cores are spaced by distance I (in units of a), 
and all doped holes are captured by the vortices, then 
x = fy, and the magnetic periodicity is 21, i.e. <5 = 




FIG. 2: (Color online) Dependence of neutron peak splitting 
parameter 8 on the doping concentration x. Filled circles 
represent the experimental results in Lai.e-iNdo^Sr^CuO^, 
empty circles — in L^-^Sr^CuOzr^. Solid line represents 
Eq.©i small-dash line Eq.© with x c o = 0.35, large-dash 
line is the linear dependence 5 = x. 



Therefore, 




This dependence is compared with the experimental re- 
sults for I^-zSr^CuO^ and Lai.g-xNdo^Sr^CuO^ 
in Fig. dJ 

It is unrealistic to expect that all doped holes are cap- 
tured inside spin-poor regions of spin superstructures. 
Some may remain delocalized in spin-rich regions. In this 
case, the above expressions would overestimate 8 both for 
stripes and for spin vortices. Therefore, the S-vs.-x scal- 
ing as such cannot conlusively confirm or rule out any of 
the two propositions. Yet, the experimental data shown 
in Fig. [5] appear to lend somewhat stronger support to 
the spin vortex scaling. 

One can further assume that only holes added in excess 
of certain concentration xq are captured by the vortices. 
In this case 



S=^-(x-x Q ). (3) 

An example of such dependence with xq — 0.035 is plot- 
ted in Fig. 2. Concentration xq itself can and should fur- 
ther depend on the net doping level x. However, it is not 
unreasonable to assume that it stays approximately the 
same in not too broad range of x, e.g. between x = 0.1 
and x = 0.15. 

Absence of higher order spin harmonics. 

Another outstanding issue that needs to be understood 
within the frameworks of both stripes and spin vortex 
interpretations is the lack of evidence for the higher or- 
der spin harmonics^. In general, the higher order peaks 
have intensity smaller than the main ones, but given the 



3 



strength and the sharpness of the latter, it is difficult to 
understand the complete lack of evidence for the former. 
Higher order harmonics have been observed for stripe- 
ordered nickelates^i. 

Spin vortex lattice does offer a good argument against 
the presence of higher order spin harmonics: As ex- 
plained below, one cannot sustain a non-collinear pattern 
of staggered spin polarizations, if only a few harmonics 
are added to the main ones. One needs a simultaneous 
contribution of many additional harmonics of comparable 
amplitude, in order to satisfy a local stability condition. 
In such a case, the harmonics with the lowest suscepti- 
bility (highest energy cost) would control the higher or- 
der spin response, and therefore, this response would be 
anomalously low. 

I assume the nearest-neighbor exchange interaction 
characterized by energy 



NN 

E = J S Sjj • S, ( 

ij,m>i,n>j 



(4) 



where (ij) and (run) are the pairs of square lattice in- 
dices, Sij are the local staggered spin polarizations, and 
J is the exchange constant. Superscript "NN" indicates 
that sites (mn) are the nearest neighbors of sites (ij)- 
Given Eq.(|4]), each spin experiences local field: 



NN 

fry = ^ . 



(5) 



A non-collinear spin configuration is stable locally, 

^ij,xhij : y Sij^yHij x 0, (6) 



when hij || Sy or 



where subscripts x and y indicate the projections of re- 
spective vectors. If this condition is violated, the local 
polarization cannot be static, because it would start pro- 
cessing around the local field direction. 

It has been shown in Ref£ that the above condition is 
satisfied for a superposition of any two non-helical sinu- 
soidal harmonics. Here I consider what happens in the 
case of an arbitrary superposition of helical or non- helical 
harmonics with in-plane spin polarization having form 



a 

= S Q (r jj ) , 



.fwcos (q Q 

S a yCOS (q Q 



fax) 
fay) 



(7) 



where a is the index of spin harmonics, q Q the wave 
vectors; S ax , S ay and fax, fay are, respectively, the 
amplitudes and the phases for the x- and y-components 
of a given harmonic; S Q (r.y) denotes the entire spin field 
of one harmonic. 

According to Eq.©, the local field corresponding to 
the static spin pattern given by Eq.© is 

hjj = -2 J ^2 [ cos (q« ' a ) + cos • b)] S Q (r.y ), (8) 

a 

where a and b are the primary translation vectors along 
the two principal directions of the underlying square lat- 
tice. [Note: If the exchange with more remote neighbors 
were included in energy ((4]), it would add terms such as, 
e.g. cos (q Q • (a + b)) or cos (q Q • 2a) to the prefactor in 
the square bracket in Eq.@. The presence of these terms 
would not change the conclusions that follow.] 

The substitution of expressions ([7]) and ([5]) into the 
local stability condition © gives 



V] [cos (qg • a) + cos (q,g • b) - cos (q a • a) - cos (q Q • b)] 

a,/3>a 

X [S ax S/3 y COS (q 

a * ^ij ~T~ (pax 

) cos (qg • Tij + cpp y ) - S ay Sp x co$ (q a ■ r y + f ay ) cos (q^ • r y + fp x )\ = 0. (9) 
I 



There are two particularly simple situations, when con- 
dition (j9]) is satisfied. The first one corresponds to all 
harmonics collinear along the same direction. In this 
case the second square bracket in each term of Eq.© 
equals zero. Both stripe and grid interpretations fulfill 
this condition, and, therefore, the local stability does not 
prevent either stripes or grid from having any combina- 
tion of higher order harmonics collinear with the main 
ones. 

The second simple situation corresponds to the first 
square bracket in Eq. © being equal to zero for all pairs of 
a and (3. This condition is satisfied, for example, for any 
two harmonics commensurate or incommensurate, helical 



or not, provided q Q and q^ are rotated with respect to 
each other by 90 degrees, and |q Q | = |qg|. Therefore, any 
non-collinear structure involving only two suchharmonics 
qi and q2 (given earlier) is locally stable. 

The wave vectors of potential higher order harmonics 
have form 

/ 2ir 2tt \ 
q a = ± — n, ± — m ) ; with n, m — 0, 1, 2, 3, 4; (10) 
\ 8a 8a J 

except for the pairs (n = ±l,m = 0) and (n = 0, m — 
±1) corresponding to the wave vectors ±q x and q2, re- 
spectively. 

When two non-collinear harmonics qi and q2 are al- 
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ready present, none of the higher order harmonics (|T0|) 
would make the first square bracket in Eq.© equal to 
zero for the term involving this higher order harmonic 
and either qi or q2. In this case, condition © cannot 
be satisfied on term-by-term basis for every pair of har- 
monics a and j3. However, it is still possible to make 
the total sum @ identically zero for every value of r^. 
The Fourier transform of 8 x 8 spin superstructure has 
maximum 64 wave vectors. When all possible harmonics 
contribute to the spin superstructure, there should exist 
infinitely many ways of fulfilling condition because 
two polarization directions of the contributing harmonics 
(x and y) imply 128 adjustable parameters (amplitudes 
and phases) to cancel 64 Fourier components of the sum 
(J9j> . The non-collinear superstructure considered in Ref. 5 
consisting of disconnected collinear magnetic clusters, is 
one such a locally stable example of adding higher or- 
der harmonics to the ones at qi and q2- There, in fact, 
every cluster can have an arbitrary polarization with re- 
spect to the others and still remain locally stable. Such 
a superstructure does contain sizable contributions from 
many higher order spin harmonics, none of which has 
been observed experimentally. 

It is now important to realize, that one cannot simply 
add a small number of higher order harmonics to the two 
at qi and q2 and still preserve the local stability. The 
reason is the following: 

Each (a,/3)-term in sum ([9]) contains a product of 



cos(q Q • Tij + tp___) and cos(q / 3 • r.^ + ip,,,) and, there- 
fore, contributes to the Fourier components of the sum 
at q a + qp and q a — qp. The contributions to a given 
Fourier component can thus come from different pairs 
(a, /?), which need to cancel each other. Therefore, there 
should be either zero or at least two such contributions 
to a given Fourier component. With too few harmon- 
ics added to qi and q2, there will necessarily be values 
of q a ±qp, where only one (a,/3)-term contributes, and 
therefore the local stability condition cannot be satisfied. 

I was not able to establish rigorously the minimum 
number of higher order harmonics, to cancel only the 
leading contributions to sum ([9]) from pairs involving qi 
or q2, but from experimenting with various possibilities, 
the set harmonics located at 

/ 2tt 2tt \ 
q a = ±—n, ±—m ; with n, m — 0, 1, 2; (11) 
\ 4a 4a / 

appears to be the most economic proposition. It is clear 
in any case, that one needs to introduce a large number 
of harmonics of comparable amplitude spanning the en- 
tire Brillouin zone, which brings us back to the argument 
made earlier that the susceptibility to such a configura- 
tion should be anomalously low. 

The author is grateful to E. Carlson, T. Egami, M. 
Fujita, H. Koizumi and G. Seibold for discussions and 
communications related to this work. 
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